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The effective potential for the composite fields responsible for chiral symmetry breaking in weakly 
coupled QED in a magnetic field is derived. The global minimum of the effective potential is found 
ON to acquire a non-vanishing expectation value of the composite fields that leads to generating the 

dynamical fermion mass by an external magnetic field. The results are compared with those for the 
Nambu-Jona-Lasinio model. 
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I. INTRODUCTION 



The problem of symmetry breaking and vacuum instability under the influence of external fields has attracted a 
lot of attention in the past few years Recently, some progress was attained in understanding the mechanism 

of the so-called catalysis of chiral symmetry breaking p|-|lO[] (for some earlier studies see also fll]-|l3|). This problem 
provides an especially interesting example of vacuum engineering: i.e., manipulating external fields to alter the 
■^j- ' symmetry properties of the vacuum. It is also one of a small number of tractable non-perturbative problems in (3+1)- 
7-H , dimensional quantum field theories, and as such may provide insight into the non-perturbative vacuum structure of 
' other theories such as quantum chromodynamics. 

Here we continue the study of this effect in weakly coupled quantum electrodynamics (QED). In particular, we derive 
the expression for the effective potential of composite fields responsible for chiral symmetry breaking by applying the 
method of jl^Jl^], originally introduced in the gauged Nambu-Jona-Lasinio (NJL) model and QED3. 
■ In this paper we are interested in constructing the effective potential for the composite local fields a(x) — 
(0|?/'(a;)V ; ( :z; )|0) and ir(x) = (0|?/'(^)*75''/'( a; )|0)- I n general, this problem requires consideration of the correspond- 
ing generating functional, 

Z{J S , J p ) = exp [iW(J s , J p )] = J Y[ (dip(x)di(j(x)dA tl (x)) exp i J d 4 x (C + J s (x)y 7 )(x)ip(x) + J p {x)tp(x)i-y 5 il)(x)) 

• <-H 1 x 

X- (i) 

where C is the Lagrangian density of massless QED in an external magnetic field. Then the calculation of the effective 
action reduces to obtaining the Legendre transform of the functional W(J S , J p ) with respect to the external sources, 

I>,tt) = W{J S , Jp) - J d 4 x [J s (x)a(x) + J p (x)n(x)} . (2) 
Here the sources are functions of fields, which are obtained by inverting the following expressions: 

, , sw 



The latter, in turn, mean that 



SJ s (x) SJp(x) 



(5r (5r 

-J s (x), —— = -J p (x). (4) 



Sa(x) Stt(x) 
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II. THE EFFECTIVE POTENTIAL FOR COMPOSITE FIELDS 



Due to the presence of chiral symmetry, the effective potential is a function of the only invariant, p = \J a 2 + 7r 2 . 
Thus, it is sufficient to keep only one of the sources to be non-zero (we choose J s = J, and J p = 0). At the end, the 
effective potential of one variable a could be promoted to the complete expression by the substitution \a\ — > p. 

Since the constant source (this is just enough for obtaining the effective potential) enters the expression for the 
generating functional in Eq. (yj) in exactly the same way as a bare mass would (toq — — J), the calculation could be 
considerably simplified |l4| , |l5| ]. 

Indeed, from the definition of the a — (0\ipip\0)j and the first relation in Eq. (^), we obtain [Tifl 

So (J) 

w(J)= J (0|V^|0),/^dS , (5) 
where, by definition, Y.q is the value of self-energy function at p = 0, and 

W(J) = w(J) { d A x. (6) 



The chiral condensate is given by the trace of the fermion propagator. In weakly coupled QED, it is adequate to use 
the so-called lowest Landau level approximation. In this approximation, the fermion propagator (with the momentum 
dependent dynamical mass £(p||)) is given by 

S^( P ) = i exp (Vx/ 2 ) C 1 ~ *rV) , (7) 



where / = 1/ -\/|e-B| is the magnetic length and, without loss of generality, we assume that sign(e-B) = 1. Here we 
have adopted the same notations as in Refs. 0,0. Also, it is reasonable to use a sharp momentum cut-off at around 



the scale of magnetic field, i.e. A ~ y|e5?l- The latter means that all the "bare" quantities in such an approach are 
defined at the characteristic scale of magnetic field. 

The main argument in support of such an approximation is the presence of the large Landau gap of order ^(ei?! 
in the energy spectrum for the modes from the higher Landau levels and the absence of any (or the presence of 
a small one after dynamical fermion mass is generated ) for the modes in the lowest Landau level. Besides, the 
infrared region with k « \feB has been shown to play a dominant role in the catalysis of chiral symmetry breaking 
by an external magnetic field . 

In the lowest Landau level approximation, from Eq. (Q) we derive 

, . — . . 1 f d(k 2 )E(k) 

- J (8) 

where we switched to Euclidean space. We would like to mention here that the above expression for the chiral 
condensate can also be obtained by making use of the E p representation (see Appendix B in the second paper of 
Ref. [El). The Schwinger-Dyson equation for the self-energy function, in the same approximation, reads [^) (compare 
also with Eq. (106) in the second paper of Ref. [Q) 

EM--J+— f^M. ?^pH^/2) 
[P) ~ + 2W k 2 + T?{k)J (k-p) 2 + q ■ W 

o 

Note that here we added the bare mass term tuq — — J, which is absent in 

To find a non-perturbative solution to the above equation for the self-energy, we use the method of Q (see Ap- 
pendix C of the second paper). As is common in solving such integral equations, we approximate it by its linearized 
version. While looking for a solution in the form £ = S(fc 2 ), we can perform the angular integration in Eq. (||) exactly. 
Therefore, we arrive at 

A 2 oo „ . 

a f d(fc 2 )S(fc 2 ) f dqexp(-ql 2 /2) 



E(p 2 ) = - J+ — / \' ^' / - v ' ' (10) 

K ' 2* J k 2 + ^ 2 J ^(4:2+^2 + g )2_ 4fc2p 2 v i 
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Finally, after approximating the kernel by its asymptotes in two regions k 2 <C p 2 and k 2 3> p 2 and dropping the 
exponential in the second integral^ (which is justified in our theory with a sharp cut-off at the scale of magnetic field), 
we arrive at the following integral equation, 



x 

where we use the dimcnsionless quantities 



M(y) , « 

(11) 



V + M 2 



M = El, j = Jl, x = p 2 l 2 , y = k 2 l 2 . (12) 



It is straightforward to check that the integral equation (11) is equivalent to the following differential equation 



d 2 M(x) 1 dM{x) a M(. 



x 



+ -— ir^ + ^z — rrF2T = o, (13) 



cPx x dx 2n x(x + M 
with a solution subject to the following (infrared and ultraviolet) boundary conditions: 



(xM(x))\ x=o =0, ( M {x)-xhx{x)^^\ 



-j- (14) 



x=l 



Note also that using the definition of the chiral condensate in Eq. (||), we arrive at the following expression: 

(15) 



/«, t 1 f , M {y) 1 f dM(x) 

0# j = — — ttt / — -fro = ~ — x— 
o 



2 = 1 



where, to be consistent with the approximations to the Schwinger-Dyson equation, we again used the linearized version 
instead of the exact relation. 

The solution to Eq. ([l3]) that is regular at the origin (in order to satisfy the infrared boundary condition in Eq. (|Ti|)) 
is the hypergeometric function Q , 

M{x)=M F(iv,-iv;l;-j^), v = ^ . (16) 

The ultraviolet boundary condition leads to the following relation between Mq and the external source: 

M F(iv,-iu;l ] -^) = -j. (17) 

Now we are in a position to derive the generatingfunctional and the effective potential. In weakly coupled QED, the 
solution for the self-energy is such that M <C 1 (jj^l- Therefore, using the asymptotic behavior of the hypergeometric 
function J16|, we obtain 



..-Mo^a.MnUeM-*.,,). (is, 

1 V 77 T 



(0\W\0)J - ^V^^sin^ln^ + ^M) , (19) 



where we have introduced the following function (which is of order v 2 in the coupling) 

T(l-2«/) 
,T 2 (1- iu) 



. Til - 2iv)\ 



1 We would like to thank V.P. Gusynin for suggesting this approximation. 
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Making use of the expressions in Eqs. (|lj) and ([l9|), we calculate the generating functional, 



M 2 0) 



w(J) = - 



tanh(7ri/) 
8tt 3 ^ 4 
Mq tanh(7w) 



d/isin (^ln(/i) + < E>(^)) [cos (fln(/x) + - 2^ sin (i^ln(/x) + 

[sin (2 i /ln(Af 2 ) + 2$(z/)) - 2i/] . 



16ttV 2 Z 4 

Then the effective potential V(p) is given by the following parametric representation: 



V(M ) = j 
p(M ) = I 



dw ( dj 



dM \dM 
dw ( dj 
ZMq \d~Mo~ 



Mq tanh(7ri/) 
16tt 3 i/ 2 Z 4 



sin 



>ln(M 2 ) +2$(i/)) +2i/] , 



47T 2 J// 3 



Mq / tanh(-7w) 



sin(j/ln(M 2 ) + <J>(t/))| 



(21) 

(22) 
(23) 



Here, as we argued above, we substituted the chiral invariant p instead of |er|. 

We notice that the effective potential given by Eqs. (p^), ( p3[ ) is a multivalued function. This, at first sight, unusual 
property of the potential originates from the presence of an infinite tower of excited resonances with the same quantum 
numbers as those of a and it particles. Among all of these branches, we have to choose the only one which corresponds 
to the stable states. 

The gap equation, derived from this effective potential, is equivalent to the equation j = 0. The latter leads to the 
infinite set of solutions 



M< n) =exp (~|-(2n + l)) , n = 0,l,2,. 

'tanh(7w) / 7r . 

— i — ^exp- — (277 + 1) 



» - 



1 



4tt 2 v1 3 



(24) 
(25) 



Note that the solutions with negative values of n are all spurious, since they correspond to values of the dynamical 
fermion mass larger than the cut-off (we remind that A = l/l here) . All the solutions in Eq. (|2^) are local minima of 
the effective potential. Indeed, 



d 2 V(p) 



dp 2 



1 dj ( dp 



I dM \dM 



8ttVZ 2 > 0. 



(26) 



However, the only solution that represents the stable composite a and 7r particles is that with 77 = Q. We could 
also check that the value of the effective potential at p^ is lower than any of the values at the other minima, as it 
should be. In figure |l| the potential is plotted in the physically interesting region in the vicinity of p(°\ In this region, 
the potential as a function of p is approximately given by the following expression: 



1 - In 



P 



,(0) 



v(p) ~ -2^vzy 

Clearly, the dynamical fermion mass m^ yn that corresponds to the global minimum p^ is 

m dyn = M^/l ~ y/\eB\ exp (~\J^j i 
which is consistent with the result obtained in Refs. 



(27) 



(28) 



III. CONCLUSION 



In this paper we derived the effective potential of the composite fields responsible for chiral symmetry breaking in 
weakly coupled QED in a constant magnetic field. The global minimum of the effective potential is obtained and is 
found to acquire a non-vanishing expectation value of the composite fields that leads to generating the dynamical 
fermion mass by an external magnetic field. Comparing the results in QED with those in the NJL model S, we 
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notice a qualitative difference in their effective potentials. While in the NJL model the potential is a single valued 
function, in QED it is multivalued. The latter property reflects the presence of a tower of excited states with the 
same quantum numbers as those of the Nambu-Goldstone bosons. From the physical point of view, this difference 
appears due to the long-range interaction in QED in contrast to the short-range interaction in NJL model. 

To calculate the full low-energy effective action, we also need the kinetic terms for these fields. Applying the 
general method of Rcfs. Jl7],[L8|], the kinetic part of the effective Lagrangian density can be expressed in terms of the 
propagators of composite scalar and pseudoscalar fields, A s (g) and A p (g): 



1 

4 



A 



dqndqx 



q=0 



1 d 2 (ApHg)) 



d^Ti(x)d\T:(x). 



(29) 



q=0 



The propagators, entering the right hand side of this formula, depend on the non-perturbative fcrmion propagator (the 
solution to the Schwinger-Dyson equation) and the scalar and pseudoscalar vertex functions p!7|^8[ . The calculation 
of these vertex functions requires the solution of an additional non-perturbative self-consistency equation analogous 
to the Schwinger-Dyson equation, which is beyond the scope of the present paper. We expect that the forms of the 
kinetic terms and the corresponding dispersion relations should be qualitatively similar to those derived in the case 
of the NJL model ||. Important points are: 1) there is no mass gap in the pseudoscalar dispersion relation, since 
the pseudoscalar particle is a Goldstone boson. 2) For the scalar particle there is a mass gap, and we expect that, as 
in the NJL model, the coefficients of the kinetic terms will be different for the directions parallel and perpendicular 
to the applied magnetic field. The dependence of the dispersion law on the transverse components is expected to be 
strongly suppressed, but it is nonetheless extremely important, playing the key role in preventing the chiral symmetry 
breaking from being washed out. 
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FIG. 1. The effective potential at v = 0.2 in the vicinity of the n — minimum. The units are chosen so that I = 1. The x 
axis represents the chiral invariant p. 
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